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Attorney Docket No: NAVI-009/02US PATENT 



DECLARATION 
As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated next to my name. 

I believe I am the original, first and sole inventor (if only one name is listed below) or an 
original, first and joint inventor (if plural names are listed below) of the subject matter which is 
claimed and for which a patent is sought on the invention entitled: 

3 METHOD FOR SCREENING AND PRODUCING COMPOUND LIBRARIES 

ffl the specification of which: 
\f\ [X] is attached hereto. 

5-5 
r = % 

hj [] was filed on , and identified as Attorney Docket No. NAVI-009/02US. 

^s, [] was filed on , as Application Serial No. 

"g [] the amendment(s) of which were filed on . 

I herebv state that I have reviewed and understand the contents of the above-identified 
specification, including the claims, as amended by any amendment referred to above. 

I acknowledge the duty to disclose information which is material to the patentability of 
this application in accordance with Title 37, Code of Federal Regulations, Section 1.56. 

I hereby claim foreign priority benefits under title 35, United States Code, Section 1 19 of 
any foreign application(s) for patent or inventor's certificate listed below and have also identified 
below any foreign application for patent or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 

Prior Foreign Application(s) (Country) (Number) (Day/Month/Year Filed) Priority Claimed (Yes/No) 



I hereby claim the benefit under Title 35, United States Code, § 119(e) of any United 
States provisional application(s) listed below. 
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60/100.290 

(Application Number) 



September 14. 1998 
(Filing Date) 



60/100.224 September 14. 1998 

(Application Number) (Filing Date) 

60/109. 2 32 November 18. 1998 

(Application Number) (Filing Date) 

60/109 ,234 Novem b er IS, 1998 

(Application Number) (Filing Date) 

I hereby claim the benefit under Title 35, United States Code, Section 120 of any United 
'States application(s) listed below and, insofar as the subject matter of each of the claims of this 
application is not disclosed in the prior United States application in the manner provided by the 
first paragraph of Title 35, United States Code, Section 112, 1 acknowledge the duty to disclose 
material information as defined in Title 37, Code of Federal Regulations, Section 1.56(a) which 
occurred between the filing date of the prior application and the national or PCT international 
filing date of this application: 

Appl. Ser, Nq, Filing Date Status fPatu/Pend./AbandV 



I direct that correspondence concerning this application be directed to 

COOLEY GOD WARD LLP 
Attention: Patent Group 
Five Palo Alto Square 
3000 El Camino Real 
Palo Alto, California 94306-2155 
Telephone (650) 843-5000. 
Facsimile (650) 857-0663 

I hereby declare that all statements made herein of my own knowledge are true and that all 
statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may jeopardize the validity of the application or any 
patent issued thereon. 
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Full name of sole or First inventor: i^Grass, George M. 

Inventor's signature Date f/l j/ff 

Residence: / Tahoe City, California 

Citizen of: United States of America 

Post Office Address: 1 506 Juniper Mountain Road, P.O. Box 1 242 

Tahoe City, California 96 1 45 

Full name of sole or Second inventor: Leesman, Glen D. 

Inventor's signature l^*7^*6£L~ Date / St / / 

Residence: Hamilton, Montana 

Citizen of: United States of America 

Post Office Address: 1 86 Nighthawk Lane 

Hamilton, Montana 59840-9307 

Full name of sole or Tl/ir^ inve#ftr: Nojfis, Daniel A. 

Inventor's signature Ib&Utl/ StQWs^ Date & Jl j^^f 

Residence: San Diego, California 

Citizen of: United States of America 

Post Office Address: 3 145 Cowley Way, #130 

San Diego, California 921 17 

Full name of sole or Fourth inventor: Sinko, Patrick J. 

Inventor's signature ( j Date ^ 1 

Residence: ^Lebanon, New Jersey 

Citizen of: United States of America 

Post Office Address: 2 Country Place 

Lebanon, New Jersey 08833 
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Full name of sole or Fifth invent 

Inventor's signature 



Residence: 
Citizen of: 
Post Office Address: 




Wehrli, John E. 



Date. 



Mountain View, California 

United States of America 

1879 Springer Unit B 

Mountain View, California 94040 
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Attorney Docket No: NAVI-009/02US 



PATENT 



Express Mail Label Number: 
Date of Deposit: 



EM 570 539 405 US 
May 26, 1999 



I hereby certify that this paper or fee is being deposited with the United States Postal Service "Express Mail 
Post Office to Addressee" service under 37 CFR 1.10 on the date indicated above and is addressed to the 
Assistant Commissioner for Patents, Washington, DC 2023 1 . 



Date: 



May 99 



By: 



j/LJH)M//2_ s*lu GA- 



YS THE UNITED STATES PATENT AND TRADEMARK OFFICE 



In re application of George M. Grass, et al. 
Serial No.: Not yet assigned 
Filed: Herewith 



Examiner: Not yet assigned 
Art Unit: Not yet assigned 



For: 



METHOD FOR SCREENING AND PRODUCING COMPOUND LIBRARIES 



Assistant Commissioner for Patents 
BOX PATENT APPLICATION 
Washington, D.C. 20231 

POWER OF ATTORNEY BY ASSIGNEE AND EXCLUSION OF INVENTORS 

UNDER 37 CFR 1.36 AND 3.71 

Sir: 

The undersigned assignee of the entire interest in the application for 
Letters Patent identified above hereby revokes all prior appointments of attorneys and 
appoints 



Nina M. Ashton 


37,273 


Marcella Lillis 


36,583 


Alexandra J. Baran 


39,101 


Tom M. Moran 


26,314 


James A. Bradburne 


38,389 


Richard L. Neeley 


30,092 


Aaron S. Brodsky 


39,920 


Craig P. Opperman 


37,078 


Shelley P. Eberle 


31,411 


Marya A. Postner 


42,085 


Richard M. Goldman 


25,585 


Gurjeev K. Sachdeva 


37,434 


Willis E. Higgins 


23,025 


William E. Winters 


42,232 


Peter R. Leal 


24,226 


Kevin J. Zimmer 


36,977 
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all of the firm of Cooley Godward LLP, to prosecute this application and to transact all 
business in the United States Patent and Trademark Office connected therewith. This 
appointment shall be to the exclusion of the above-identified inventor(s) and any 
attorney(s) appointed by such inventor(s), in accordance with the provisions of 37 C.F.R. 
1.36 and 3.71. 

Assignee's rights are evidenced by an assignment 

[X] a copy of which is enclosed herewith. 

[] previously recorded on at reel , frame(s) . 



Please direct all telephone calls and correspondence to: 



Cooley Godward LLP 
Attention: Patent Group 
Five Palo Alto Square 
3000 El Camino Real 
Palo Alto, CA 94306-2155 
Telephone: 650-843-5000 
Facsimile: 650-857-0663 




Title: 



Senior Director. Legal Affai rs/ Corporate Secretary 



Address: 9880 Campus Point Drive. San Diego. California 92121 



Date: 



May 14, 1999 
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Attorney Docket No: NAVI-009/02US PATENT 



Express Mail Label Number: EM 570 539 405 US 

Date of Deposit: May 26, 1 999 

I hereby certify that this paper or fee is being deposited with the United States Postal Service "Express Mail 
Post Office to Addressee" service under 37 CFR 1.10 on the date indicated above and is addressed to the 
Assistant Commissioner for Patents, Washington, DC 20231. 



/Vh# 99 



Date: 

t* & irusf 7 7 By: 




IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re application of George M. Grass, et al. 

Serial No.: Not yet assigned Examiner: Not yet assigned 

Filed: May 26, 1999 Art Unit: Not yet assigned 



For: METHOD FOR SCREENING AND PRODUCING COMPOUND LIBRARIES 



Assistant Commissioner for Patents 
Washington, D.C 20231 



STATEMENT UNDER 37 CFR 3.73(b) 
ESTABLISHING RIGHT OF ASSIGNEE TO TAKE ACTION 

1 . The assignee(s) of the entire right, title and interest hereby seek(s) to take action 
in the PTO in this matter. 

IDENTIFICATION OF ASSIGNEE 

Name of Assignee: Navicyte. In corporated ; 

Type of Assignee: Corporation 



1 
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PERSON AUTHORIZED TO SIGN 

Name of Person Authorized to Sign: James A. Bradburne, Ph.D. 

Title of Person Authorized to Sign: Agent of Record 

[X] I, the person signing below, state that I am empowered to sign this statement on 
behalf of the assignee. 

BASIS OF ASSIGNEE'S INTEREST 

Ownership by the assignee is established as follows. A chain of title from the inventor(s) 
to the current assignee is shown below: 

1. From: George M. Grass. Glen D. Leesman. Daniel A. Norris, 

Patrick J. Sinko and John E. Wehrli 

To: Navicvte, Incorporated __ 



Recordation Date: 12/22/98 Reel: 9673 Frame: 0774-0777 



[ ] Copies of the documents in the chain of title are attached. 



Cooley Godward LLP Respectfully submitted, 

Attn: Patent Group Cooley Godward llp 

Five Palo Alto Square 
3000 El Camino Real 

Palo Alto, CA 94306-2155 n C\s<Ai jlJ?< 

Tel: (650) 843-5000 ^ ^HMr^ 
ew r*.sn\ zen-ntti James A. Bradb 




Fax: (650) 857-0663 James A - Bradbunie, Ph.D. 

' Reg. No. 38,389 

JAB:pkp 
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PATENT APPLICATION 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
In re the Application of: 

GRASS et al . Art Unit: 1 631 

Application No.: 09/320,069 Examiner: M. Moran 

Filed: May 26, 1999 Attorney Dkt. No.: 109904-00027 

For: METHOD FOR SCREENING AND PRODUCING COMPOUND LIBRARIES 

REVOCATION OF POWER OF ATTORNEY AND NEW APPOINTMENT 

Assistant Commissioner for Patents 

Washington, D.C. 20231 November 5, 2001 

Sir: 

Dr. Christian Kilger, the undersigned, certifies that Lion Bioscience AG is 
the assignee of the entire right, title and interest in U.S. Patent Application Serial No. 
09/320,069, filed May 26, 1999, by virtue of an assignment from NaviCyte, Inc., a copy 
of which is attached hereto. The assignment to NaviCyte was recorded in the Patent 
and Trademark Office at Reel 9673, Frame 0774-0777 on December 22, 1998. 

Further, Lion Bioscience AG, as assignee of the entire interest in and to 
the above-identified United States patent application hereby revokes all powers of 
attorney previously given and appoints Arent Fox Kintner Plotkin & Kahn, 1050 
Connecticut Avenue, Suite 600, Washington, DC, 20036-5339, a firm composed of: 



TECH/48828. 1 



Charles M. Marmelstein, Reg. No. 25,895; Robert B. Murray, 
Reg. No. 22,980; George E. Oram, Jr., Reg. No. 27,931; 
Douglas H. Goldhush, Reg. No. 33,125; Richard J. Berman, 
Reg. No. 39,107; Murat Ozgu, Reg. No. 44,275; Robert K. 
Carpenter, Reg. No. 34,794; Gregory B. Kang, Reg. No. 
45,273; Rustan Hill, Reg. No. 37,351; Kevin F. Turner, Reg. 
No. 43,437; Rhonda C. Barton, Reg. No. P47.271 and Hans 
J. Crosby, Reg. No. 44,634, Brian A. Tollefson, Reg. No. 
46,338, Lynn D. Anderson, Reg. No. 46,412, David D. Dzara, 
Reg. No. 47,543; Laurence J. Edson, Reg. No. 44,666; 
Michael A. Steinberg, Reg. No. 43,160; and Lynn A. Bristol, 
Reg. No. 48,898. 



as principal attorneys to prosecute said application and to transact all business in the 
Patent and Trademark Office connected therewith. 



The undersigned has reviewed all of the appropriate documents and, to the 



best of the undersigned's knowledge and belief, title is in the assignee identified above. 



The undersigned (whose title is supplied below) is empowered to sign this 



By the provisions of 28 U.S.C. §1746, I hereby declare under penalty of 



perjury under the laws of the United States of America that the foregoing is true and 



paper on behalf of the assignee. 



correct. 



Executed on: November 5, 2001 

(date) 



Name: 



Dr. Christian Kilger 



TECH/48796. 1 



Title: 



Signature: 




ATTACHED SHEET 



U.S. Patent 5,183,760 

U.S. Patent 5,591,636 

U.S. Serial Number 09/320,069 

U.S. Serial Number 09/320,270 

U.S. Serial Number 09/320,545 

U.S. Serial Number 60/224,106 



U.S. Patent 5,599,688 

U.S. Patent 6,146,883 

U.S. Serial Number 09/320,372 

U.S. Serial Number 09/320,371 

U.S. Serial Number 09/320,544 



ASSIGNMENT 



WHEREAS, NaviCyte, Inc. (hereinafter ASSIGNOR) a company having its principal 
office and place of business at Sparks, Nevada, USA owns the inventions, and the patents 
and patent applications in various countries directed thereto, as set forth on the 2 

page attachment hereto. 

AND, WHEREAS, LION Bioscience AG (hereinafter ASSIGNEE), a company having 
its principal office and place of business at Heidelberg, Germany, is desirous of acquiring 
all interest therein. 

NOW, THEREFORE, in consideration of One Dollar ($1.00) and other good and 
valuable consideration, the receipt and sufficiency of which are hereby acknowledged, the 
above said ASSIGNOR by these presents does sell, assign and transfer unto the said 
ASSIGNEE, its successors, assigns and legal representatives, the full and exclusive right, 
title and interest in and to the said inventions in the United States and other countries 
including the ful and extensive right, title and interest world-wide in and to the patents and 
patent applications on the attachment hereto. 



NaviCyte, Inc. 





Witness: 



APPENDICIES 



Appendix 1: Abbreviation Key for Mass- Volume Model 



Abbreviation 



Kf sd = associated rate constant for stomach and duodenum 



Ka dj = associated rate constant for duodenum and jejunum 



Ka ji = associated rate constant for jejunum and ileum 



Ka ie = associated rate constant for ileum and colon 



Ka co = associated rate constant for colon and excretion 



SD trans = transfer rate between stomach and duodenum 



DJ trans = transfer rate between duodenum and jejunum 



JL trans = transfer rate between jejunum and ileum 



IC trans = transfer rate between ileum and colon 



Waste = transfer rate between colon and excretion 



pH s = pH stomach 



pH s2 = pH duodenum 



pH s3 = pH jejunum 



pH s4 = pH ileum 



pH s5 = pH colon 



sol profile = solubility profile for stomach 
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sol profile 2 = solubility profile for duodenum 



sol profile 3 = solubility profile for jejunum 



sol profile 4 = solubility profile for ileum 



sol profile 5 = solubility profile for colon 



stom ka = associated rate constant for stomach compartments 1 and 2 



duo ka = associated rate constant for duodenum compartments 1 and 2 



Jej ka = associated rate constant for jejunum compartments 1 and 2 



II ka = associated rate constant for ileum compartments 1 and 2 



Colon ka = associated rate constant for colon compartments 1 and 2 



SA stom = surface area of stomach 



SA duo = surface area of duodenum 



SA jej = surface area of jejunum 



SA il = surface area of ileum 



S A colon = surface area of colon 



Perm stom = permeability of stomach 



Perm duo = permeability of duodenum 



Perm jej = permeability of jejunum 



Perm il = permeability of ileum 



Perm colon = permeability of colon 
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Ka sd = associated rate construct for stomach fluid absorption 



Ka du = associated rate construct for duodeunm fluid absorption 



Ka je = associated rate construct for jejunm fluid absorption 



Ka il = associated rate construct for ileunm fluid absorption 



Ka co = associated rate construct for colon fluid absorption 



Note: other abbreviations adhere to above descriptors and are self 
explanatory 
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Appendix 2: Equations, Parameters and Values For Mass-Volume Model 

amt_plasma(t) = amt_plasma(t - dt) + (trans_2 1 + ka - elimination - trans_12) * dt 
INIT amt_plasma = 0 

INFLOWS: 

trans_2 1 = k21*comp_2 
ka = tot_abs_rate 

OUTFLOWS: 

elimination = amt_plasma*k_elim 
trans_12 = kl2*amt_plasma 

blood_side_col(t) = blood_side_col(t - dt) + (colon_ka_5) * dt 
INIT blood_side_col = 0 

INFLOWS: 

colon_ka_5 = IF Vol_colon*sol_profile_5 >=Colon THEN Colon*SA_colon*perm_colon*3600 
ELSE Vol_colon*sol_profile_5*SA_colon*perm_colon*3600 
blood_side_dou(t) = blood_side_dou(t - dt) + (duo_ka) * dt 
INIT blood_side_dou = 0 

INFLOWS: 

duo_ka = IF Vol_duod*sol_profile_2 >= duodenum THEN 
duodenum*SA_duo*perm_duo*3600ELSE Vol_duod*sol_profile_2*SA_duo*perm_duo*3600 
blood_side_il(t) = blood_side_il(t - dt) + (Il_ka) * dt 
INIT blood_side_il = 0 

INFLOWS: 

Il_ka = IF Vol_ileum*sol_profile_4 >=Ileum THEN Ileum*SA_Il*perm_Il*3600 ELSE 
Vol_ileum*sol_profile_4*SA_Il*perm_Il*3600 
blood_sideJej(t) = blood_sideJej(t - dt) + (Jej_ka) * dt 
INIT blood_sidejej = 0 

INFLOWS: 

Jej_ka = IF VolJej*sol_profile_3 >=Jejunum THEN Jejunum*SAjej*permJej *3600 ELSE 
Vol Jej *sol_profile_3 * S A Jej *perm Jej * 3600 
blood_side_sto(t) = blood_side_sto(t - dt) + (stom_ka) * dt 
INIT blood_side_sto - 0 

INFLOWS: 

stom_ka = IF Vol_stom*sol_profile >= Stomach THEN Stomach*SA_stom*perm_stom*3600 
ELSE Vol_stom*soI_profile* SA_stom*perm_stom*3600 
Colon(t) = Colon(t - dt) + (ICjrans - Waste - colon_ka_5) * dt 
INIT Colon = 0 
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INFLOWS: 

ICtrans = ka_ic*Ileum 

OLTFLOWS: 

Waste = ka_col*Colon 

colon_ka_5 = IF Vol_colon*sol_profile_5 >=Colon THEN Colon* SA_colon*perm_colon* 3600 
ELSE Vol_colon* sol_profile_5 * S A_colon*perm_colon* 3 600 
comp_2(t) = comp_2(t - dt) + (trans_12 - trans_21) * dt 
INIT comp_2 = 0 

INFLOWS: 

trans_12 = kl2*amt_plasma 

OUTFLOWS: 

trans_21 = k21*comp_2 

duodenum(t) = duodenum(t - dt) + (SD_trans - duojca - DJ_trans) * dt 
q INIT duodenum = 0 

yQ INFLOWS: 

CO SDjrans = if Stomach >0 then kf_sd* Stomach else 0 

Ul OUTFLOWS: 

UJ duo_ka = IF Vol_duod*sol_profile_2 >= duodenum THEN 
duodenum*SA_duo*perm_duo*3600 ELSE Vol_duod*sol_profile_2*SA_duo*perm_duo*3600 

f D J_trans = ka_dj * duodenum 

£ excretion(t) = excretion(t - dt) + (vol_cw) * dt 

?T INIT excretion = 0 

'fl INFLOWS: 
£7 vol_cw = Vol_colon*ka_col 
'~ excretion_2(t) = excretion_2(t - dt) + (Waste) * dt 
INIT excretion_2 = 0 

INFLOWS: 

Waste = ka_col*Colon 

Ileum(t) = Ileum(t - dt) + (JLjrans - IC_trans - Il_ka) * dt 
INIT Ileum = 0 

INFLOWS: 

JL_trans = kaJi*Jejunum 

OUTFLOWS: 
IC_trans = ka_ic*Ileum 

Il_ka = IF Vol_ileum*sol_profile_4 >=Ileum THEN Ileum*SA_Il*perm_Il*3600 ELSE 

Vol_ileum*sol_profile_4*SA_Il*perm_Il*3600 

Jejunum(t) = Jejunum(t - dt) + (DJ_trans - JL_trans - Jej_ka) * dt 
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INIT Jejunum = 0 
INFLOWS: 

DJtrans = kadj* duodenum 

OUTFLOWS: 

JLtrans = kaJi*Jejunum 

Jejka = IF VolJej*sol_profile_3 >=Jejunum THEN Jejunum*SAJej*permJej 

Vol Jej *sol_profile_3 *SA Jej*permjej *3600 

serosal col(t) = serosal_col(t - dt) + (Adsorp_col - col_secretion) * dt 

INIT serosalcol = 0 

INFLOWS: 

Adsorp_col = PULSE(1.67,0,.l)+0*Vol_colon*ka_co 

OUTFLOWS: 
col_secretion = 0 

serosal_dou(t) = serosal_dou(t - dt) + (Adsorp_Duo - duo_secretion) * dt 
INIT serosal_dou = 0 

INFLOWS: 

Adsorp_Duo = PULSE(1 0.82,0,. l)+0*Vol_duod*ka_du 
OUTFLOWS: 

duo_secretion = PULSE(1 0.82,0,. 1) 

serosaMll(t) = serosal_ill(t - dt) + (Adsorpt_ill - ile_secretion) * dt 
INIT serosal_ill = 0 

INFLOWS: 

Adsorpt_ill = PULSE(8.83,0,.10)+0*Vol_ileum*ka_il 
OUTFLOWS: 

ile_secretion = PULSE(1.50,0,.l) 

serosaljej(t) = serosal Jej(t - dt) + (Adsorpjej - jej_secretion) * dt 
INIT serosaljej = 0 

INFLOWS: 

Adsorpjej = PULSE(1 5.76,0,. l)+0*VolJej*kaJe 
OUTFLOWS: 

jej_secretion = PULSE(2.67,0,. 1 ) 

serosaI_sto(t) = serosal_sto(t - dt) + (Adsorp_Stom - Stom_Secretion) * dt 
INIT serosal_sto = 0 

INFLOWS: 

Adsorp_Stom = 0*Vol_stom*ka_sd 
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OUTFLOWS: 

StomSecretion = PULSE(1 6.67,0,. 1 ) 

Stomach(t) = Stomach(t - dt) + (- SDjrans - stom_ka) * dt 

INIT Stomach = 1000 

OUTFLOWS: 

SDtrans = if Stomach >0 then kf_sd*Stomach else 0 

stom_ka = IF Vol_stom*sol_profile >= Stomach THEN Stomach*SA_stom*perm_stom*3600 
ELSE Volstom * sol_profile * S Astom * perm_stom * 3 600 
totaldrugabsorbed(t) = total_drug_absorbed(t - dt) + (tot_abs_rate) * dt 
INIT total_drug_absorbed = 0 

INFLOWS: 

totabsrate = stom_ka+duo_ka+Jej_ka+Il_ka+colon_ka_5 
Total_Elimination(t) = Total_Elimination(t - dt) + (elimination) * dt 
INIT Total_Elimination = 0 

INFLOWS: 

elimination = amt_plasma*k_elim 

Vol_colon(t) = Vol_colon(t - dt) + (vol_ij + col secretion - vol cw - Adsorp col) * dt 
INIT Vol_colon = 0 

INFLOWS: 

volij = Vol_ileum*ka_ic 
colsecretion = 0 

OUTFLOWS: 

vol_cw = Vol_colon*ka_col 

Adsorp_col = PULSE(1.67,0,.l)+0*Vol_colon*ka_co 

Vol_duod(t) = Vol_duod(t - dt) + (vol_sd + duo secretion - voil dj - Adsorp Duo) * dt 
INIT Vol_duod = 0 " ~ 

INFLOWS: 

vol_sd = kf_sd*Vol_stom 
duosecretion = PULSE(1 0.82,0,. 1) 

OUTFLOWS: 

voil_dj = Vol_duod*ka_dj 

Adsorp_Duo = PULSE(10.82,0,.l)+0*Vol_duod*ka_du 

Vol_ileum(t) = Vol_ileum(t - dt) + (volji + ile_secretion - Adsorpt ill - vol ij) * dt 
INIT VoMleum = 0 " 

INFLOWS: 

volji = VolJej*kaJi 

ile_secretion = PULSE( 1.50,0,.!) 
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OUTFLOWS: 

Adsorpt_ill = PULSE(8.83,0,.10)+0*Vol_ileum*ka_il 
volij = Vol_ileum*ka_ic 

Voljej(t) = VolJej(t - dt) + (voil_dj + jej_secretion - vol Ji - Adsorp jei) * dt 
INITVolJej=0 

INFLOWS: 

voildj = Vol_duod*ka_dj 
jej_secretion = PULSE(2.67,0,.l) 

OUTFLOWS: 
volji = VolJej*kaJi 

Adsorpjej = PULSE(l5.76,0,.l)+0*VolJej*kaJe 

Vol stom(t) = Vol_stom(t - dt) + (Stom_Secretion - vol_sd - Adsorp Stom) * dt 
INIT Vol_stom = PULSE(8.33,0,-1) 

INFLOWS: 

Stom_Secretion = PULSE(1 6.67,0,. 1) 

OUTFLOWS: 

vol_sd = kf_sd*Vol_stom 

Adsorp_Stom = 0*Vol_stom*ka_sd 

conc_plasma = (amt_plasma/volume)*mg_to_ug 

kl2 = .839 

k21=.67 

kaco = 1 

kacol = 3 

ka_dj = 3 

kadu =1 

kaic = 3 

ka_il = 8.83 

kaje = 1 

kaji = 3 

kasd = 1 

kf_sd = 2.8 

kelim = .161 

mg_to_ug= 1000 

permcolon = 3.80e-6 

permduo =1.1 0e-6 

permll = 4.06e-6 

permjej =2.17e-6 

permstom = 1.10e-6 

ph_s=1.5 

ph_s_2 = 6.6 

ph_s_3 = 6.6 
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ph_s_4 = 7.5 
ph_s_5 = 6.6 
SA_colon= 138 
SA_duo= 125 
SA_I1 = 102 
SAJej = 182 
SA_stom = 50 
volume = 4* 19200 
sol_profile = GRAPH(ph_s) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65) (4 50 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00 3.65)' 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_2 = GRAPH(ph_s_2) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4 50 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3 65)' 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_3 = GRAPH(ph_s_3) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4 50 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65) 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_4 = GRAPH(ph_s_4) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00. 3.65), (4 50 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65X (8.00, 3.65)' 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_5 = GRAPH(ph_s_5) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4 50 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65)' 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
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Appendix 3: Abbreviation Key For GI Model 



The legend/key has been divided into sub-sections corresponding to the sub-sections of the 
model diagram. 

Numbered suffixes (1, 2, 3, 4, 5, 6) have been assigned to distinguish between intestinal regions 
(stomach, duodenum, jejunum, ileum, colon, and waste, respectively). 

1 - stomach 

2 - duodenum 

3 - jejunum 

4 - ileum 

5 - colon 

6 - waste 

For example, VOL 1 is the volume in the stomach, MASS 3 is the insoluble mass in the jejunum. 
In the equations, COMP 1 indicates the stomach, COMP 2 the duodenum, COMP 3, the jejunum, 
etc. 

Ghosts are listed under the sub-section containing the original reservoir, flow regulator, or 
converter. 

Abbreviations listed in italics are regionally dependent and set up as arrays to allow independent 
values for each intestinal region. 

In general, ADJ as a prefix indicates a calculated parameter value (ADJ - adjusted), while ADJ 
as a suffix indicates an adjustment parameter (ADJ = adjustment). 



Intestinal model 

Reservoirs/Compartments 

Fluid volume absorbed 
Fluid volume 

Mass of drug contained with a formulation or controlled release 
device 

Insoluble mass of drug (not contained within the formulation or 
controlled release device) 
Soluble mass of drug 
Mass of drug absorbed 

Flow regulators 



VOL ABS 

VOL 

CREL 

MASS 

SOL 

ABSORPTION 
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REABS Rate of water absorption 

VOL OUT Fluid volume transit rate 

CR OUT Formulation or controlled release device transit rate 

CR INPUT Drug release rate from formulation or controlled release device 

MASS OUT Insoluble drug mass transit rate 

DISSPRECIP Dissolution rate 

SOL OUT Soluble drug mass transit rate 

FLUX Absorption rate 



ADJ PARMS (Adjustment Parameters) 

VOL ADJ Fluid volume absorption adjustment parameter 

DISS ADJ Dissolution rate adjustment parameter 

TRANSIT ADJ Transit time adjustment parameter 

SA ADJ Surface area adjustment parameter 

FLUX ADJ Passive Absorption adjustment parameter 

EFFLUX ADJ Efflux or secretion adjustment parameter 

CARRIER ADJ Active absorption adjustment parameter 



PARMS (Parameters) 



VOL PARM 
SURFACE AREA 
DOSE 

INIT VOLUME 
TIME IN HOURS 
pH 

PARACELLULAR 



Fluid volume absorption rate constant 
Surface area available for absorption 
The administered dose of drug 
The administered volume of water or fluid 
A clock 

The physiological pH value 

A user controlled switch used to adjust absorption based on 
absorption mechanism 



TRANSIT TIME 



TRANSFERS 
CUMUTT 

ADJ TRANSIT TIME 
USER TT INPUT 



GI transit rate constant 
Cumulative transit time 

Adjusted GI transit time incorporating adjustment parameter and 
user input 

User controlled adjustments to the GI transit time 



OUTPUT CALCULATIONS 

ABSORBED TOTAL Total mass of drug absorbed (sum of ABSORPTION 1 ...5) 



107. 



Fraction or the dose absorbed into portal vein x 1 00 

Total absorption rate (sum of FLUX 1 ...5) 

Cumulative drug mass dissolved 

Cumulative drug mass released from formulation 

Sum of DISS PRECIP 1...5 

Summ ofCR INPUT 1...5 



PERMEABILITY CALCULATION 

Adjusted permeability ncorporating all transport mechanisms and 
relevant adjustment parameters 
Active or carrier-mediated absorptive permeability 
Constant from the Michaelis-Mentin type permeability equation for 
active transport 

Passive permeability after regional correlation calculation (same as 
PASS PE if regional correlation is not used) 
Passive permeability entered by user 

A logical function used in determining the regional correlation * 
A logical function used in determining the regional correlation 

SOLUBILITY CALCULATION 

USER pH User supplied pH value for which a solubility value is available 

USER SOLUB User supplied solubility value corresponding to the USER pH value 

ADJ SOLUB Solubility calculated (if necessary) at the appropriate pH value 

using the entered USER pH and USER SOLUB values 



CONTROLLED RELEASE CALCULATION 

CR RATE The instantaneous release rate from the formulation 

CR DOSE The total dose contained with the formulation 

CR AT TIME The cumulative drug mass release profile 

CR AT LAST The cumulative drug mass release profile 

Note: CR AT TIME holds the value at the current time value (t), CR AT LAST holds the value 
at the immediately proceeding time value (t-1) 



CONC CALCULATION 

CONCENTRATIONS The dissolved drug concentration 



FDp% 

FLUX TOTAL 
CUM DISS 
CR Release 
CUM DISS RATE 
CR cumrate 



ADJ PERM 

ACT PE 
Km 

REGIONAL 

PASS PE 
RC 

RCSUM 
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DISSOLUTION CALCULATION 

PRECIP Precipitation rate constant 

DISSOL Dissolution rate constant 

ADJ DISS PRECIP Adjusted rate constant incorporating PRECIP, DISSOL and 

calculated concentration 
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Appendix 4: Equations, Parameters and Values For GI Model 



□I' ADJ PARMS 

G CARRIER. ADJ[COMPS] = 0 
G DISS_ADJ[COMP_1] = 1 
G DISS_ADJ[COMP~2] = 1 
O DISS_ADJ[COMP_3] = 1 
3 DISS_ADJ[COMP~4] = 1 
G DISS_ADJ[COMP_5] = 1 
Q EFFLUX_ADJ[COMPS] = 1 
O FLUX_ADJ[COMP_1] = 1 
G FLUX_ADJ[COMP_2] = 10 
G FLUX_ADJ[COMP_3] = 8 
3 FLUX_ADJ[COMP_4] = 2 
G FLUX_ADJ[COMP_5] = 1 
G SA_ADJ[COMP_1] = 1 
O SA_ADJ[COMP_2] = 1 
G SA_ADJ[COMP_3] = 1 
G SA_ADJ[COMP_4] = 1 
G SA_ADJ[COMP_5] = 1 
G TRANSIT_ADJ[COMP_1] = 1 
G TRANSIT_ADJ[COMP_2] = 1 
Q TRANSIT_ADJ[COMPl3] = 1 
O TRANSIT_ADJ[COMP_4] = 1 
G TRANSIT_ADJ[COMP_5] = 1 
G VOL_ADJ[COMP_1] =~1 
O VOL_ADJ[COMP_2] = 1 
O VOL_ADJ[COMP_3] = 1 
G VOL_ADJ[COMP_4] = 1 
G VOL_ADJ[COMP_5] = 1 
(g CONC CALCULATION 

G CONCENTRATIONS[COMP_1] = if VOL_1=0.0 then 0 else if 

ADJ_SOLUB[COMP_1]<SOL_1/VOL_1 then ADJ_SOLUB[COMP_1] else SOL 1/VOL 1 + 

0*(SOL_2+SOL_5+SOL_3+SOL_4+VOL_3+VOL_2+VOL_4+VOL 5) 
G CONCENTRATIONSlCOMP_2] = if VOL_2 = 0.0 then 0 else if (VOL_2<1e-6 AND SOL_2<1e-7) 

then 0 else if ADJ_SOLUB[COMP_2]<SOL_2/VOL_2 then ADJ SOLUB[COMP 21 else 

SOL_2A/OL_2 

+0*(SOL_1 +SOL_5+SOL_3*SOL_4+VOL_3+VOL 1 +VOL_5+VOL_4) 
O CONCENTRATIONS[COMP_3J = if VOL_3 - 0.0 then 0 else if (VOL_3<1e-6 AND SOL_3<1e-7) 
then 0 else if ADJ_SOLUB[COMP 3]<SOL 3/VOL 3 then ADJ_SOLUB[COMP 3] else 
SOL_3A/OL_3 

♦0*(SOL_1 +SOL_2+SOL_4+SOL_5+VOL_5+VOL_4+VOL_1 +VOL 2) 
G CONCENTRATIONS[COMP_4] ="if VOL_4 = 0.0 then 0 else if (VOL_4<1e-6 AND SOL 4<le-7) 
then 0 else if ADJ_SOLUB[COMP 4]<SOL_4/VOL 4 then ADJ SOLUBfCOMP 4] else" 
SOL_4/VOL_4 

+0*(SOL_1 +SOL_2+SOL_3+SOL_5+VOL_1 +VOL_2+VOL_3+VOL_5) 



I/O 



Z C0NCENTRATI0NS[C0MP_5] s if V0L_5 * 0.0 then 0 else if (VOL 5<1e-6ANDSOL 5<le 7> 
SOL 5^0L 5 J - S ° LUBlC0MP - 5 ^ 

+0*(SOL_H-SOL_4+SOL_3+SOL_2+VOL 3+VOL 1+VOL 2+VOL 41 
,gi CONTROL RELEASE CALCULATION " " " ~ 

Q CR_DOSE = 0 

C CR_RATE = (CR_AT_TIME-CR AT_LASTr20*(CR_DOSE7100) 
O CR_AT_LAST = GRAPH(TIME-DT) 

(0 00, 0.00), (0.25. 17.7), (0.5. 31.5), (0.75. 42.2). (1.00, 50.6). (1.25. 57.1) (1 50 62 1) (1 75 
66.1), (2.00. 69.2). (2.25. 71.6). (2.50. 73.4). (2.75. 74.9), (3.00 76.0) (3 25 76 9W3 50 7 b\ 
(3 75, 78.1). (4.00. 78.5). (4.25. 78.9). (4.50. 79.1). (4.75. 79.3) (5 00 79 5)' (5 25 79 6)' (S 50 
79.7). (5.75. 79.7). (6.00. 79.8). (6.25. 79.8), (6.50. 79.9). (6.75 79 9)' (7 00 79 7» n f 5 
(7-50. 80.0). (7.75. 80.0). (8.00. 80.0). (8.25. 80.0). (8.50. 80.0) (8 75 900) (9 00 80 0) (9 25 
M°,°i 1™°' 80 0) ' (9 75 ' 80 0) ' (10 0 ' 80 0) ' < 10 - 3 ' 80 °)- < 105 ' 80.0) (10 8.' 80 0) (11 0 80 0 ' 

□ 50 0) { ' 8 °- 0) ' (11 - 8, 80 - 0)l (12 - 0, 80 0)> (12 - 3 ' 8 °- 0) - {12 - 5 ' 80 -°>- < 12 8 80 OK fl 3 0 

'5 0 CR_AT_TIME = GRAPH(TIME) 

*j (0.00. 0.00). (0.25. 17.7). (0.5. 31.5). (0.75. 42.2). (1.00. 50.6). (1.25. 57.1) (1 50 62 1) (1 75 

~ « ^ £° 0, 69 2)> (2 - 25, 71 - 6) ' (2 ' 50 ' 734 >' < 2 - 75 ' ™>- (3-00 76.0). (3.25 76 9) (3 50 77 6 ) 
m % 7 *> 78 , 1 J. ( 4 -0°. 78.5). (4.25. 78.9). (4.50. 79.1). (4.75. 79.3). (5.00 79 5) (5 25 79 6) (5 50 
h ^ (6.00. 79.8). (6.25. 79.8). (6.50. 79.9). (6.75, 79.9) (7.00 79 0) (7 25 79 9) 

Si m ffS' 7 ft 75 ' 800) ' (8 °°' 800) ' (8 ' 25 ' 80 °)' ( 850 ' 8 °0). (8.75. 80.0) (9 00 80 0) (9 25 
80.0), (9.50. 80.0), (9.75. 80.0). (10.0, 80.0). (10.3. 80.0). (10.5. 80.0) (10 8 80 0) (11 0 80 0) 
U (1 o 1 o 3. 80.0). (11.5. 80.0). (11.8. 80.0). (12.0. 80.0). (12.3. 80.0). (12 5. 80 0).' (12 ^ 80 0)! (13 A 

jy) DISSOLUTION CALCULATION 

Z O ADJ_DISS_PRECIP[COMP_1]aifVOL 1*0 then 0 else if 
^ (SOL_1 A/OL_1 <ADJ_SOLUB[COMP_lf) then 

f (pjSf^OMP-^'D'SS.ADJICOMP.IJ-MASS.I^ADJ.SOLUBfCOMP 1]-SOL 1/VOL 1))else 
((SOL.1/VOL_1).ADJ_SOLUB[COMP_1]rPRECIP[COMP 1]+ ~ " VUL - 1J,e,Se 

0*(MASS_UMASS_2+MASS_3+MASS_4+MASS 5+SOL~1+SOL 2+SOL 3+SOL 4+SOL 5-V 
OL_1+VOL_2+VOL_3+VOL 4+VOL 5) . - - °^ L - J OUL -« &UL_b+V 

O ADJ_DISS_PRECIP[C0MPJ2] = if VOL 2=0 then 0 else if 
(SOL_2/VOL_2<ADJ_SOLUB[COMP 2)) then 

;PJ 8S °4 C ° MP - 2 1* DIS ^^^ 2J-SOL 2/VOL 2)) else 

((SOL_2A/OL_2)-AOJ_SOLUB[COMP 2]rPRECIP[COMP_2] ' ^L.^eise 

+0'(MASS_UMASS_2+MASS_3+MASS_4+MASS 5+SOL 1+SOL 2+SOL 3+SOL 4+SOL 5+V 
OL_1+VOL_2+VOL_3*VOL 4+VOL 5) ----- v 

O ADJ_DiSS_PRECIP[COMP3] = if VOL 3*0 then 0 else if 
(SOL_3/VOL_3<ADJ SOLUB[COMP 3]) then 

(P^f 0 ^^- 31 *^^^ 3]-SOL 3/VOL 3)) else 

((SOL.3A/OL_3)-AOJ_SOLUB[COMP 3]rPREc7P[COMP 3] " " ,,e ' Se 

+0*(MASS_1+MASS_2+MASS_3+MASS_4+MASS 5*SOL HSOL 2+SOL 3+SOL 4+SOL 5+V 
OL_1+VOL_2+VOL_3+VOL_4+VOL_5) " " " ~ ~ " 



Hi 



ADJ_DISS_PRECIP[C0MP_4] s if VOL 4=0 then 0 else if 
i SOL_4/VOL_4<ADJ_SOLUB[COMP_4f) then 

(DlSSOL[COMP_4]-DISS_ADJ[COMP_4]*MASS_4-(ADJ SOLU8[COMP 41-SOL 4/VOL 4)) else 
((SOL_4/VOL_4)-ADJ_SOLUB[COMP 4])'PRECIP[COMP 4] " ~ 

♦0-(MASS_UMASS_2+MASS_3+MASS_4+MASS_5+SOL 1+SOL 2+SOL 3+SOL 4+SOL 5W 
OL_1 +VOL_2+VOL_3+VOL_4+VOL_5) " " " ~ 

ADJ_DISS_PRECIP[COMP_5] * if VOL_5=0 then 0 else if 
(SOL_5/VOL_5<ADJ_SOLUB[COMP_5])then 

(DISSOL[COMP_5rDISS.ADJ[COMP_5]*MASS 5*(ADJ SOLUBfCOMP 51-SOL 5/VOL 5)) else 
((SOL_5/VOL_5)-ADJ_SOLUB[COMP_5]rPREcTP[COMP_5) " ~ " 

+0'(MASS_1+MASS_2+MASS_3+MASS_4*MASS_5+SOL 1+SOL 2+SOL 3+SOL 4+SOL 5+V 
OL_1+VOL_2+VOL_3+VOL_4+VOL_5) " 





DISSOL(COMP_1] 


= 1 


^— - 


DISSOL{COMPl2] 


= 1 




DISSOL[COMP_3] 


= 1 




DlSSOL[COMP_4] 


= 1 


5 


DISSOL[COMP~5] 


= 1 


i" 


PRECIP[C0MPJ] 


= 10 


0 


PRECIP[C0MP_2] 


= 10 


0 


PRECIP[C0MP_3] 


= 10 


<■ — 


PRECIP[COMP_4] 


= 10 


c 


PRECIP[COMP~5] 


= 10 


INPUTS 




INTESTINAL MODEL 





□ ABSORPTIONS (t)* ABSORPTION 1 (t - dt) + (FLUX 1) • dt 
IN IT ABSORPTIONS =0 

INFLOWS: 

FLUX_1 = 

CONCENTRATIONStCOMP_irADJ_PERM[COMP_irSURFACE_AREA[COMP_1] 

□ ABSORPTION_2(t) * ABSORPTION_2(t - dt) ♦ (FLUX.2) * dt 
INIT ABSORPTION 2*0 

INFLOWS: 

FLUX_2 = 

CONCENTRATIONS[COMP_2] , ADJ_PERM[COMP_2] # SURFACE_AREA[COMP_2] 

□ ABSORPTION_3(t) * ABSORPTION 3(t - dt) + (FLUX.3) * dt 
INIT ABSORPTION 3*0 

INFLOWS: 
■£> FLUX_3 = 

CONCENTRATIONStCOMP_3]*ADJ_PERM[COMP_3]'SURFACE_AREA[COMP_3] 

□ ABSORPTION_4(t) - ABSORPTION_4(t - dt) ♦ (FLUX_4) * dt 
INIT ABSORPTION 4*0 

INFLOWS: 



FLUX_4 = 

CONCENTRATIONS[COMP_4rADJ_PERM[COMP_4]«SURFACE_AREA[COMP_4] 

d ABS0RPTI0N_5(t) = ABSORPTION 5(t - dt) + (FLUX 5) • dt 
I NIT ABSORPTION_5 = 0 
INFLOWS: 

FLUX_5 = iftime<32then 

CONCENTRATIONS[COMP_5]*ADJ_PERM[COMP_5]*SURFACE AREA[COMP 5H32-ti 
me)/48*(VOL_5/1 7.2) else 0 ~ ~ 

U2 C_REL_1 (t) = C_REL_1 (t - dt) + (- CR_OUT_1 - CR INPUT 1 ) * dt 
INIT C_REL_1 = CR_DOSE 
OUTFLOWS: 

CR_OUT_1 = IF TIME >= CUMU_TT[COMP 1] THEN C REL_1 '10000 ELSE 0 
CR_INPUT_1 =jfTIME>CUMU_TT[COMP ijthen 0 else~CR_RATE 
C C_REL_2(t) = C_REL_2(t-dt) + (CR_OUT 1-CR OUT 2 - CR INPUT 2)*dt 
INIT C_REL_2 = 0 ----_/ 

INFLOWS: 

CR_OUT_1 = IF TIME >= CUMU_TT[COMP_1] THEN C REL 1*10000 ELSE 0 
OUTFLOWS: " " 

-ff CR_OUT_2 = IF TIME >= CUMU_TT[COMP 2] THEN C REL_2*10000 ELSE 0 

CR_INPUT_2 = ifTIME>CUMU_TT[COMP_2Tthen0else~CR_RATE 

□ C_REL_3(t) = C_REL_3(t-dt) + (CR_OUT 2-CR OUT 3-CR INPUT 3) * dt 
INIT C_REL_3 = 0 " " ~ 

INFLOWS: 

CR_OUT_2 = IF TIME >= CUMUJT[COMP 2] THEN C REL 2*10000 ELSE 0 
OUTFLOWS: " ~ 

# CR_OUT_3 = IF TIME >= CUMU_TT[COMP 3] THEN C_REL 3*10000 ELSE 0 
CRJNPUT_3 = if TIME > CUMU_TT[COMP 3] then 0 else CR^RATE 

□ C_REL_4(t) = C_REL_4(t - dt) ♦ (CR_OUT 3 - CR OUT 4 - CR INPUT 4) ' dt 
INIT C_REL_4 = 0 

INFLOWS: 

CR_OUT_3 = IF TIME >= CUMU TT[COMP 3] THEN C REL 3*10000 ELSE 0 
OUTFLOWS: ~ ~ 

•§> CR_OUT_4 = IF TIME >= CUMU_TT[COMP_4] THEN C REL_4*1 0000 ELSE 0 
CRJNPUT4 = if TIME>CUMU_TT[COMP 4? then 0 else CR_RATE 

□ C_REL_5(t) = C_REL_5(t.dt) + (CR_6uT 4-CR~OUT 5-CR INPUT 5) * dt 
INIT C_REL 5 = 0 ~ " 

INFLOWS:" 

CR_OUT_4 = IF TIME >= CUMU TT[COMP_4] THEN C_REL 4*10000 ELSE 0 
OUTFLOWS: 

CR_OUT_5 = IF TIME >= CUMU_TT[COMP_5] THEN C REL 5*10000 ELSE 0 
^ CR_INPUT_5 = ifTIME>CUMU TT[COMP 5] then 0 else CR RATE 

□ C_REL_6(t) = C_REL_6(t - dt) ♦ (CR_OUT_5) * dt" 
INIT C_REL 6 = 0 

INFLOWS:" 

CR_OUT_5 = IF TIME >= CUMU_TT[COMP_5] THEN C_REL_5*10000 ELSE 0 



Il3 



CZ MASS_1(t) = MASS_1(t-dt) + (CRJNPUT 1 - MASS OUT 1 - DISS PRECIP 1) 
INIT MASS_1 = DOSE " 
INFLOWS: 

CR_INPUT_1 = if TIME>CUMU_TT{COMP_1] then 0 else CR RATE 
OUTFLOWS: 

«jf MASS.OUTJ =MASS_1*TRANSFERS[COMP 1] 

# DISS_PRECIP_1 = ADJ_DISS_PREClP[COMPJ] 

□ MASS_2(t) = MASS_2(t - dt) ♦ (MASS.OUT 1 ♦ CRJNPUT 2 - MASS OUT 2 - 
DISS_PRECIP_2)*dt " " 

INIT MASS_2 = 0 
INFLOWS: 

MASS_OUT_1 =MASS_1TRANSFERS[COMP 1] 
CR_INPUT_2 * if TIME>CUMU_TT[COMP_2] then 0 else CR RATE 
OUTFLOWS: 

MASS_OUT_2 = MASS_2*TRANSFERS[COMP 2] 
DISS_PRECIP_2 = ADJ_DISS_PRECIP[COMP3] 

□ MASS_3(t) = MASS_3(t - dt) + (CRJNPUT 3 + MASS OUT_2 - MASS OUT 3 - 
DISS_PRECIP_3) * dt ~ " 

INIT MASS 3 = 0 
INFLOWSP 

CR_INPUT_3 = ifTIME>CUMU_TTlCOMP 3] then 0 else CR_RATE 
MASS_OUT_2 = MASS_2 # TRANSFERSICOMP 2] 
OUTFLOWS: 

MASS_OUT_3 = MASS_3TRANSFERS[COMP 3] 
DISS_PRECIP_3 = AOJ_DISS_PRECIP[COMP~3] 

□ MASS_4(t) = MASS_4(t-dt) + (CR INPUT_4 ♦ MASS_OLIT 3 - MASS OUT 4- 
DISS_PRECIP 4) • dt ~ 

INIT MASS_4 = 0 
INFLOWS: 

•$ CRJNPUT_4 = ifTIME>CUMUJT[COMP 4] then 0 else CR_RATE 
MASS_OUT_3 = MASS_3TRANSFERS[COMP 3] 
OUTFLOWS: 

# MASS.OUT4 = MASS_4TRANSFERS[COMP_4] 
DISS_PRECIP_4 - ADJ_DISS_PRECIP[COMP_4] 

□ MASS_5(t)*MASS_5(t-dt) + (CRJNPUT 5 ♦ MASS OUT 4 -MASS OUT 5- 
DISS_PRECIP_5) * dt ~ ~ ~ 

INIT MASS.5 = 0 
INFLOWS: 

CR_INPUT_5 * if TIME>CUMU_TT[COMP_5] then 0 else CR RATE 
MASS_OUT_4 = MASS 4TRANSFERS[COMP_4] 

OUTFLOWS: 

^ MASS_OUT_5 » if time>4 then MASS_5*TRANSFERS[COMP 5] else 0 

# DISS_PRECIP_5 * ADJ_DISS_PREC7P[C0MP 5] 

□ MASS_6(t) * MASS 6(t • dt) ♦ (MASS OUT 5) • dt 
INIT MASS.6 = 0 

INFLOWS: 



MASS.0UT.5 = if tim >4 th n MASS_5*TRANSFERS[COMP_5] else 0 
— S0L_1(t) = S0L_1(t - dt) + (DISS.PRECIP 1 - SOL OUT 1 - FLUX 1 r dt 
INIT SOL_1 =0 - 
INFLOWS: 

DISS_PRECIP_1 *ADJ_DISS_PRECIP[COMP 11 
OUTFLOWS: 

^ SOL_OUT_1 = SOL_TTRANSFERS[COMP 1] 
FLUX_1 = 

CONCENTRATIONS[COMP_1]«ADJ_PERM[COMP_1]'SURFACE_AREA[COMP_1] 

O SOL_2(t) = SOL_2(t - dt) + (SOL_OUT_1 + DISS.PRECIP 2 • SOL OUT 2 - FLUX 2) * dt 
INIT SOL_2 = 0 ~ 

INFLOWS: 

SOL_OUT_1 aSOL_1*TRANSFERS[COMP 1] 
DISS_PRECIP_2 « ADJ DISS_PRECIP(COMP_2] 
OUTFLOWS: 

SOL_OUT_2 = SOL_2TRANSFERS[COMP_2] 

# FLUX_2 = 

CONCENTRATIONS[COMP_2rADJ_PERM[COMP_2]*SURFACE_AREA[COMP_2] 

□ SOL_3(t)*SOL_3(t-dt) + (DISS_PRECIP_3 + SOL OUT 2 -SOL OUT 3 -FLUX 3) • dt 
INIT SOL_3 = 0 " ~ - ; 

INFLOWS: 

DISS_PRECIP_3 = ADJ_DISS_PRECIP[COMP_3] 
SOL_OUT_2 = SOL_2TRANSFERS[COMP_2] 

OUTFLOWS: 

^ SOL_OUT_3 = SOL_3TRANSFERS[COMP_3] 
FLUX_3 = 

CONCENTRATIpNS(COMP_3]*ADJ_PERM[COMP_3]*SURFACE_AREA[COMP_3] 

□ SOL_4(t)*SOL_4(t-dt) + (DISS_PRECIP 4 ♦ SOL OUT 3 -SOL OUT 4 -FLUX 4) • dt 
INIT SOL_4 = 0 - > 

INFLOWS: 

DISS_PRECIP_4 = ADJ_DISS_PRECIP[COMP_4] 

# SOL_OUT_3 = SOL 3*TRANSFERS[COMP 3] 
OUTFLOWS: 

«§> SOL_OUT_4 a SOL 4TRANSFERS[COMP 4] 
FLUX_4 = 

CONCENTRATIONS(COMP_4] , ADJ_PERM[COMP_4] # SURFACE_AREA[COMP_4] 

C SOL_5(t) = SOL_5(t - dt) + (DISS_PRECIP 5 ♦ SOL OUT 4 - SOL OUT 5 - FLUX 5) • dt 
INIT SOL_5 = 0 ~ " 

INFLOWS: 



US 



ft DISS_PRECIP_5 = ADJ_DISS_PRECIP(C0MP 5] 
^ S0L_0UT_4 = SOL_4T~RANSFERS[COMP_4]~ 

OUTFLOWS: 

°& SOL_OUT_5 = if time>4 then SOL_5*TRANSFERS[COMP_5] else 0 
ft FLUX_5siftime<32then 

CONCENTRATIONS[COMP_5]*ADJ_PERM[COMP 5]'SURFACE AREAfCOMP 51'(32-ti 
me)/48*(VOL_5/1 7.2) else 0 " ~ 

SOL_6(t) = SOL_6(t - dt) ♦ (SOL OUT_5) # dt 

INIT SOL_6 = 0 
INFLOWS: 

ft SOL_OUT_5 = if time>4 then SOL_5TRANSFERS[COMP_5] else 0 
VOL_1(t) = VOL_1(t ■ dt) ♦ (- REABS.1 - VOL.OUT 1 ) • dt 
INIT VOL_1 = INIT_VOLUME 
OUTFLOWS: 
ft REABS_1 s VOL_1*VOL_PARM[COMP_1] 
VOL_OUT_1 = VOL_1*TRANSFERS[COMP 1] 
VOL_2(t) * VOL_2(t - dt) ♦ (VOL_OUT_1 - VOL_OUT_2 - REABS 2) * dt 
INIT VOL 2 = 0 



INFLOWS: 

X ft VOL_OUT_1 =VOL_rTRANSFERS[COMP_1] 

J OUTFLOWS: 

W VOL_OUT_2 = VOL_2*TRANSFERS[COMP 2] 

hj «#> REABS.2 = VOL_2'VOL_PARM[COMP_2] " 

y □ VOL_3(t) = VOL_3(t - dt) ♦ (VOL.OUT 2 - VOL OUT 3 - REABS 3) * dt 

« . INIT VOL_3 = 0 ~ " 

U. INFLOWS: 

U ft VOL_OUT_2 = VOL_2TRANSFERS[COMP 21 

fli OUTFLOWS: 

£ # VOL_OUT_3 * VOL_3TRANSFERS[COMP_3] 

^ ft REABS.3 = VOL_3*VOL_PARM[COMP 3] 

VOL_4(t) * VOL_4(t - dt) + (VOL.OUT 3 - VOL~OUT 4 . REABS 4) • dt 

INIT VOL_4 = 0 " ~ 

INFLOWS: 

It VOL_OUT_3 = VOL_3TRANSFERS[COMP 3] 
OUTFLOWS: 

ft VOL_OUT_4*VOL_4TRANSFERS[COMP_4] 
REABS.4 = VOL_4*VOL_PARM[COMP 4] 
□ VOL_5(t) - VOL_5(t - dt) ♦ (VOL_OUT 4 - VOL~OUT 5 - REABS 5) • dt 
INIT VOL_5 = 0 ~ 
INFLOWS: 

■*> VOL_OUT_4 = VOL 4TRANSFERS[COMP_4] 
OUTFLOWS: 

VOL_OUT_5 = VOL_5TRANSFERS[COMP_5] 
ft REABS.5 = VOL_5*VOL_PARM[COMP_5J 
C VOL_6(t) = VOL_6(t - dt) ♦ (VOL OUT 5) * dt 
INIT VOL 6 = 0 



HI* 



INFLOWS: 

VOL_OUT_5 = VOL_5'TRANSFERSfCOMP_5] 

— VOL_ABS_1(t) = VOL_ABS_1(t-dt) + (REABS 1 ) • dt 
I NIT VOL_ABS_1 =0 

INFLOWS: 

REABS. 1 =VOL_rVOL_PARM[COMP_1] 
CZ V0L_ABS_2(t) = V0L_ABS_2(t - dt) ♦ (REABS 2) * dt 
I NIT V0L_ABS_2 * 0 
INFLOWS: 

^> REABS.2 = V0L_2*V0L_PARM[C0MP_2] 

□ V0L_ABS_3(t) = V0L_ABS_3(t - dt) ♦ (REABS.3) • dt 
INIT V0L_ABS_3 = 0 

INFLOWS: 

REABS_3 = V0L_3*V0L_PARM[C0MP 3] 

□ V0L_ABS_4(t) = V0L_ABS_4(t • dt) ♦ (REABS ~4) • dt 
INIT VOL ABS_4 = 0 

INFLOWS: 

REABS_4 = VOL_4*VOL_PARM[COMP_4] 
C2 V0L_ABS_5(t) = V0L_ABS_5(t - dt) ♦ (REABS 5) * dt 
INIT VOL ABS_5 = 0 
INFLOWS: 

REABS_5 = V0L_5*V0L PARMfCOMP 5] 
) MULTI DOSE CALCULATION 
| OUTPUT CALCULATIONS 

O CR_Release(t) = CR_Release(t - dt) + (CR cumrate) * dt 
INIT CR_Release = 0 
INFLOWS: 

CR_cumrate = CR_INPUT_1+CR_INPUT_2+CR_INPUT_3+CR_INPUT_4+CR_INPUT_5 

□ CUM_DISS(t) = CUM_DISS(t.dt) + (CUMM DISS RATE) * dt 
INITCUM_DISS*0 

INFLOWS: 

CUMM_DISS_RATE = 

DISS_PRECIP_1*DISS_PRECIP_2+DISS_PRECIP_3+DISS_PRECIP.4*DISS_PRECIP 

5 — _ 

O ABSORBEDJTOTAL * ABSORPTION_2MBSORPTION_3+ABSORPTION_4+ABSORPTION_5 

O PDp%aABSORBED_TOTAL/DOSE # 100 

0 FLUX_TOTAL « FLUX_2*FLUX 3+FLUX_4+FLUX_5 
PARMS 

O DOSE = 1000 

O INIT_VOLUME = 100 

O PARACELLULAR = 1 

O PH[C0MP_1] = 1.5 

O pH[COMP_2] = 5 

Q pH[COMP_3] = 6.5 
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2 pH[C0MP_4] = 7 
~ pH[C0MP_5] = 6.5 

5 toZ A ll^ [ pT- ] = " PARACELLULAR =° th6n ^SA.ADJfCOMP.l, e,e 
° 7fsZ%i^pT- 2] " " PARACELLULAR= ° then 1 50-SA_ADJ [C OMP_ 2] e,e 
C toZ A lh£™fT- 3] " " PARACELLULAR =° 1°00-SA_A DJ[ COMP_3] e<se 
G 50-8^0^^ = " PARACELLULAR=0 then 1000-SA_ADJ[CO M P_4] e.se 

3 ll^M^S^t- 51 = " PARACELLULAR= ° then 850-SA_ADJ[COMP_5] else 
O TIMEJN* HOURS = TIME 

O VOL_PARM[COMP_1] = 0*VOL_ADJ[COMP_1] 

_ C VOL_PARM[COMP_2J ' 0*VOL_ADJ[COMP 2] 

«J O VOL_PARM[COMP_3] = 1.75*VOL_ADJ[COMP_3] 

'I 0 VOL_PARM[COMP_4] = 1.75*VOL ADJ[COMP_4] 

% O VOL_PARM[COMP_5] = 0.10*VOlIaDJ[COMP_5] 
Jg] PERMEABILITY CALCULATION 
in O ACT_PE[COMPSJ = [0 , 

* = 0, 

ni o , 
r o, 

u 0] 

y, O ADJ_PERM[COMP_1] = 

ry ^i* EFFLUX -ADJ[COMP_1]))«REGIONAL[COMP 1]*FLUX ADJfCOMP 11'360<WCARRIFR 
g W|COMPjrACT.I W MPjr3aO^ 

C O ADJ_PERM[COMP_2J a 

^il E ^ UX - ADJ[C0MP - 2 J )) * REG,0N ^l C0MP 2TFLUX ADJfCOMP 21'3600+fCARRlPR 
WOMP^ACT.PEfCOMP^neOO^ 

O ADJ_PERM[COMP_3] = 

^^ EFFLUX - ADJ[C0MP - 3])),REGI0NAL I C0MP 3 J*FLUX ADJfCOMP 3]'3600+(CARRIER 
DJ[COMP_3rACT_PE[COMP_3r3600/(1 + (CONCENTRAT.ONS(COMP^(^ 

O ADJ_PERM[COMP_4] = 

{2/(1+EFFLUX_ADj"[COMP_4]))*REGIONAL[COMP 4]*FLUX ADJfCOMP 4]*3600+f CARRIER 
DJ[COMP_4rACT_PE[COMP_4r3600/(U 

O ADJ_PERM[COMP_5J = 

DJ[COMP_5] ACT_PE[COMP_5]-3600/(1+(CONCENTRATIONS[COMP 5J/(Km[COMP 5])))) " 
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r KmfCOMPS] * [1 , 

1 . 
1 , 
1 . 
1] 

" . PASS_PE[COMPS] = [0 . 
1.10E-06 . 
2.17E-06, 
4.06E-06 , 
3.80E-06 ] 

O RC[COMP_1] = PASS_PE[COMP_1J*0 
C RC[COMP_2] = IF PASS_PE[COMP_2]>0 THEN 1 ELSE 0 
O RC[COMP_3] = IF PASS_PE[COMP_3]>0 THEN 2 ELSE 0 
O RC[COMP_4] = IF PASS_PE[COMP 4]>0 THEN 4 ELSE 0 
C RC[COMP_5] = PASS_PE[COMP_5r0 
O RCSUM = RC[COMP_2]+RC[COMP_3]+RC[COMP_4] 
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_ O REGIONAL[COMP_1] = PASS PEfCOMP 1]+RCSUM*0 
a C REG"ONAL{COMP_2]«ifRCSUM=2then~ 

(EXP( -9.01 1926 + 2.594378 "LOGNd/PASS PE[COMP 21)^.065515 
*LOGN(1/PASS_PE[COMP 2]r2M-1) else 
if RCSUM=4 then 

if RCSUM=6 then 

0.5*(EXP( -9.011926 * 2.594378 'LOGNd/PASS PEfCOMP 31)-0 065515 
•LOGN(1/PASS_PE[COMP 3]r2)H-1) " -^J^ woois 

+0.5*(EXP( -21.009845 + 4.544238 *LOGN(1/PASS PEfCOMP 41) -0 140815 
•LOGN(1/PASS_PE[COMP 4])»2))*(-1 ) else " lUUM - " -° 140815 
PASS_PE[COMP_2] 
O REGIONAL[COMP_3] = ifRCSUM=1 then 

(EXP( -3.238469 ♦ 1.509131 *LOGN(1/PASS PEfCOMP 21) -0 022109 
*LOGN(1/PASS_PE[COMP 2])*2))*(-1)else~ 
if RCSUM=4 then 

£*!^ 093739 * LOGN ( 1 ^ASS^ PEfCOMP 4])»2-0 00 

23631'LOGN(1/PASSJ»E[COMP 4])*3))*(-1 ) else c^umk_4j; 2. o.OO 

if RCSUM=5 then ~ 

0.5*(EXP( -3.238469 ♦ 1.509131 *LOGN(1/PASS PEfCOMP 21) -0 022109 
*LOGN(1/PASSJ>E[COMP_2])*2))*(-1) " 

♦0.5*(EXP( -15.415683 ♦ 3.543563 *LOGN(1/PASS PEfCOMP 41) -0 100318 
*LOGN(1/PASS_PE[COMP 4])*2))*(-1) else " " 
PASS PE[COMP 3] 
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r REGI0NAL[C0MP_4] = if RCSUM=1 then 

(EXP( 14.455255 -1.264630 -L0GN(1/PASS PEfCOMP 2]) + 0 082015 
•LOGN(1/PASS_PE[COMP_2]r2)r(-1)else" 5 
if RCSUM=2 then 

(EXP( 1 1 .480333 -0.791 109 *LOGN(1/PASS_PE[COMP 3]) ♦ 0.066063 
•LOGN(1/PASS_PE[COMP_3]r2)r(-1)else " 
if RCSUM=3 then 

0.5*(EXP( 14.455255-1 .264630 *LOGN(1/PASS PE[COMP 2]) + 0 082015 
*LOGN(1/PASS_PE[COMP 2])*2))*(-1) " u.uo^uis 

+0£(EXP( 1 1 .480333 -0.791 109 *LOGN(1/PASS PE[COMP 31) ♦ 0 066063 
*LOGN(1/PASS_PE[COMP_3])*2))*(-1)else " u.uot>063 
PASS_PE[COMP_4] 
Q REGIONAL[COMP_5] = PASS_PE[COMP 5] +RCSUM*0 
©) SOLUBILIY CALCULATION 

C ADJ.SOLUB[COMPJ] = if USER_pH[COMP.1]> S pH[COMP.1] then USER.SOLUBfCOMP 1] 

1 0 e (s » 

((USER.SOLUB[COMP.3)-USER.SOLUB(COMP 2jy(USER_pHtCOMP 31-USER nHrrnuo ,1 

0 USER.SOUUB.COMP.3, 

((USER_SOLUB[COMP_4]-USER_SOLUB[COMP 3]y(USER_pH[COMP 41-USER nHfrnMP -*i 
•(PH[COMP^^^ 

SS; s h e ^ 

C elsS USER.SOLUBtCOMP., 

pH[COMP_4]-USER_pH[COMP_4])*USER SOLUB[COMP 4] else ~ 
((USER_SOLUB[COMP_4]-USER_SOLUB[COMP 3]y(USER JpHfCOMP 41-USER oHrcOMP 71 
))*(pH[COMP.4]-USER_pH[COMP 3])+USER SO~LUB[COMP 3] USER_pH[COMP_3] 

C ^^f^P-^^USER^HICOMPljrpHlCOMP sjthen USER SOLUB[COMP 31 
elserfUSER_pH[COMP_3J< P H[COMP 3] then " " LUB ^ UMK - J J 

PHfCOMP 3] 

)) (pH[COMP_3]-USER_pHfCOMP_3])+USER SOLUB[COMP 3] else ~ 

!!vn5rr7^ PH[COMP 2] 

)) (pH[COMP_3]-USER_pH[COMP 2])*USER SOLUBfCOMP 21 



USER_pH[COMPS] = [1.5 
5. 

6.5. 

7 , 

7.5] 



USER_SOLUB[COMPS] = [31 , 

3.65 , 
3.65 . 
3.65 . 
3.65] 
I TRANSIT TIME 

C ADJ_TRANSIT_TIME[C0MP_1] = .5*TRANSIT_ADJ[C0MP_1]*USER_TT_INPUT 
O ADJ_TRANSIT_TIME[COMP_2] = .25*TRANSIT_ADJ[COMP_2]*USER TT_INPUT 
O ADJ_TRANSIT_TIME[COMP_3] = 1.5TRANSIT_ADJ[COMP_3ruSERlTT_INPUT 
C ADJ_TRANSIT_TIME[COMP_4] = 1.5TRANSIT_ADJ[COMP 4]*USER TTJNPUT 
O ADJ_TRANSIT_TIME[COMP_5] = 24TRANSIT_ADJ[COMP ~5]*USER TTJNPUT 
Q CUMU_TT[COMP_1] = ADJ_TRANSIT_TIME[COMP_1] 

C CUMU_TT[C0MP_2] = ADJ_TRANSIT>IME[C0MPJ]+ADJ TRANSIT_TIME[COMP 2] 
Q CUMU_TT[COMP_3] = 

ADJ_TRANSIT_TIME[COMP_1]MDJ_TRANSIT_TIME[COMP_2]+ADJ_TRANSIT_TIME[COMP 

O CUMU_TT[COMP_4] = 

ADJ_TRANSIT_TIMEIC0MP_1]+ADJ_TRANSIT TIME[COMP 2]+ADJ TRANSIT TIMEFCOMP 

3]+ADJ_TRANSIT_TIME[COMP_4] " 
O CUMU_TT[COMP_5] = 

ADJ_TRANSIT_TIME[COMP_1]+ADJ_TRANSIT_TIME[COMP 2]+ADJ_TRANSIT TIME[COMP 

3]+ADJ_TRANSIT_TIME[COMP_4]+ADJ TRANSIT_TIME[COMP 5] 
O TRANSFERS[COMP_1] = LOGN(10yADJ_TRANSIT_TIME[COMP_1] 
O TRANSFERS{COMP_2] * LOGN(10)/ADJ_TRANSIT_TIME[COMP 2J 
O TRANSFERS[COMP_3] * LOGN(10)/ADJ_TRANSIT_TIME[COMP~3] 
O TRANSFERS[COMP_4] * LOGN(10)/ADJ_TRANSIT_TIME[COMP~4] 
O TRANSFERS[COMP_5] = LOGN(10)/ADJ TRANSIT_TIME[COMPl5] 
Q USER_TT_INPUT = 1 



